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THEORY  OF  ELECTROMAGNETIC  INSTABILITY  OF  AN  INTENSE  BEAM 
IN  A  QUADRUPOLE  FOCUSING  SYSTEM 


I.  Introduction 

Discrete  quadrupole  focusing  systems,  also  called  FODO  lattices,  have  been  used  to 
transport  charged  particle  beams  for  a  variety  of  applications.  A  FODO  lattice  is  an  alter¬ 
native  to  helical  quadrupole  windings  for  strong  focusing.  Helical  quadrupole  (stellarator) 
focusing  systems  have  been  found  to  be  subject  to  an  electromagnetic  instability,  which  we 
referred  to  as  the  three- wave  instability. 1-5  This  has  been  observed  experimentally  for  a 
325  A,  950  keV  beam  with  an  80  nsec  pulse  length.4  In  that  case,  the  stellarator  gradient 
was  400  G/cm,  the  stellarator  pitch  length  was  18  cm,  and  the  axial  magnetic  field  was  1.2 
kG.  The  growth  of  the  instability  and  subsequent  beam  loss  was  accompanied  by  3  GHz 
radiation,  as  predicted  by  theory. 

Here,  we  consider  the  corresponding  instability  for  transverse  perturbations  of  a  beam 
centroid  interacting  with  a  FODO  lattice  field  and  a  TE\\  waveguide  mode.  The  advan¬ 
tages  of  using  FODO  lattices  are:  (i)  variations  in  periodicity  are  easily  introduced,  (ii) 
variations  in  periodicity  can  change  the  character  of  electromagnetic  instability  and  re¬ 
duce  the  overall  growth  rate,  (iii)  lower  growth  rates  in  some  regimes  and  (iv)  preliminary 
particle  simulations  show  saturation  at  low  values.  Even  though  the  FODO  lattice  lacks 
a  stable  regime,  the  small  growth  rates  may  allow  a  cure  by  any  one  of  several  methods, 
such  as  loss  in  the  waveguide. 

The  dispersion  relation  may  be  solved  in  the  usual  way,  with  the  growth  rates  being 
the  imaginary  parts  of  the  roots  of  the  determinant  of  a  dispersion  matrix.  Due  to  the 
periodicity  of  the  FODO  lattice,  however,  the  dispersion  matrix  is  of  infinite  dimension. 
This  instability5,6  was  recently  analyzed  via  Floquet  theory5  and  via  approximate  disper¬ 
sion  relations  based  on  two  different  approaches.5’8  In  this  paper,  we  derive  the  dispersion 
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relation  and  show  that  the  exact  growth  rate,  in  agreement  with  the  results  from  Flo- 
quet  theory,  can  be  obtained  from  a  truncated  dispersion  matrix.  Furthermore,  we  are 
able  to  approximately  factor  a  zeroth-order  form  of  the  dispersion  relation  to  find  all  the 
beam  modes.  Identifying  the  beam  mode  that  causes  the  instability,  we  are  able  to  derive 
analytic  expressions  for  the  growth  rates  and  establish  various  regions  of  instability. 
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II.  Dispersion  Relation 


We  consider  an  alternating  gradient  quadrupole  field  (Bqx,Bqy),  where 


Bqx  =  —Bqkqf(z)y,  Bqy  =  -Bqkqf(z)x, 


(la  -  6) 


Bqkq  is  the  peak  quadrupole  field,  f(z)  is  periodic,  kq  =  2ir/\q  and  A,  is  the  period  of  the 
quadrupole  field.  The  representation  for  the  quadrupole  field  in  Eqs.  (la-b)  is  valid  near 
the  z-axis,  i.e.,  (*2  +  y 2)1/2  <<  Xq/2i r. 

In  equilibrium,  the  electron  beam  travels  along  the  axis  of  a  perfectly  conducting 
circular  waveguide  of  radius  rg  at  velocity  va  and  is  monoenergetic  with  7„  =  (1  —  /?2)-1/2, 
where  (30  =  va/c.  Both  the  beam  radius  and  beam  centroid  displacement  are  assumed  to 
be  small  in  comparison  to  the  waveguide  radius.  We  include  induced  image  charges  and 
currents  due  to  the  displaced  beam.  The  induced  electric  and  magnetic  fields  near  the 
z-axis  are 


_  nm„c‘  v  .  .  . 

—  2  —  "I" 


m0c *  V 


Bind  —  2  -T0o(ycex  —  xces), 


e  r 


(2a) 

(2b) 


where  xc(zOJz)  =  xc(za,t  =  (z-z0)/v0)  andyc(z0,z)  =  yc(z0,t  =  (z-z0)/v0)  are  the  beam 
centroid  coordinates,  v  —  u>^r2/4c2  is  Budker’s  parameter,  (i/  ~  \Ie\/170  for  an  electron 
beam  where  Ie  is  the  electron  beam  current  in  units  of  kilo-Amperes),  w*  =  4ne2ne/m0, 
ne  is  the  electron  density,  m0  is  the  electron  mass  e  is  the  elementary  charge  (assumed 
positive),  rfc  is  the  beam  radius  and  z0  is  the  axial  position  of  the  electron  at  t  =  0. 

We  expect  the  TEn  mode  to  have  the  largest  growth  rate  because  its  electric  field 
peaks  on  axis.  It  vector  potential  can  be  written  as 


A  =  [Ax(r,0,z)ex  +  Ay(r,9,z)ey)  e  xwt  -f  c.c. 
=  [A+(r,0,z)  +  A_(r,  9,  z)]  e~lu,t  +  c.c., 


(3) 


where 

A ±(r,9,z)  =  J  dk±b±n(k±)  [J0(/*nr)e±  +  J2(Hur)e±l20 e^}  exk±z, 
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A-t  are  complex  amplitudes  associated  with  the  (+)  right-handed  circularly  polarized 
(RHCP)  wave  and  (  — )  left-handed  circularly  polarized  (LHCP)  wave,  «7n  is  the  n*/,  order 
Bessel  function,  5±n  are  complex  constants,  u  is  the  radian  frequency,  k±  are  the  wave 
numbers  associated  with  the  RHCP  and  LHCP  waves,  e±  =  (e*  ±iey)/2  and  c.c.  denotes 
the  complex  conjugate.  The  boundary  condition,  requiring  the  tangential  component  of  the 
electric  field  to  vanish  on  the  perfectly  conducting  waveguide  surface,  gives  the  condition 
J[(r  —  fin Tg)  =  0,  where  Hi\Tg  is  the  the  largest  positive  zeroth  of  Bessel  function  J[  and 
1  denotes  d/dr. 

The  wave  equation  for  A  is  given  by 


c2  dt2  ) 


A  = 


(4) 


where  V^_  =  V2  —  —  is  the  transverse  Laplacian,  and  Jc  =  —  eneirr2 vc  is  the  macroscopic 
transverse  current  associated  with  the  beam  centroid.  Assuming  that  the  particle  motion 
is  driven  by  field  on  axis,  the  current  can  be  written  as 


Jc  =  -v 


m„c‘ 


6(x)S(y)( 


dxe  „ 
dt  *s 


(5) 


where  [xe(z,t),yc(z,t)]  are  the  coordinates  of  beam  centroid. 

The  wave  equation  can  be  rewritten  as 

L  [(ox  ±  iay)e~XLJt  +  c.c]  =  -4 n^6(x)6(y)^(xc  ±  iyc), 

where  ax<y  =  ( e/m0c2)AXiy  and  L  =  V2±  +  d2 /dz2  +  u;2/c2. 

The  vector  potential  in  (6)  has  the  form 

(ax  Hay)e~xu>t  =  —z  /  dkTb^1i(k^)Jo(niir)el('k:r 
m0cc  J 

m0c'  J 

After  L  operated  on  the  norm?-i/ed  vector  potential,  Eq.  (6)  can  be  written  as 

J dk;eik’*  (^-kl-^l/j  y„(m,r)6T„(*T) 


—  t  Ult 


m0c * 


(6) 


(7) 
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+  c.c. 


+  J  dk+e**-  -  4  -  )  /I(y„r)e±‘M4±11(*±) 


=  4ir-<(*)%)^(xc  ±  iye). 
c  at 


(8) 


There  are  two  equations  associated  with  (8).  We  multiply  the  equation  with  the  upper 
sign  by  f2*  d0/2ir  J*a  J0(fiur)rdr  and  the  equation  with  the  lower  sign  by  J ‘2*  ex20dO/2ir 
J^a  J2{nnr)rdr,  then  combine  the  results  to  pick  out  6_n.  Similarly,  we  multiply  the 
equation  with  the  upper  sign  by  J**  e~x20d8/2ir  /J*  J2(nnr)rdr  and  the  equation  with 
the  lower  sign  by  J*n  d8/2ir  J*a  J0(fj.nr)rdr,  then  combine  the  results  to  pick  out  6+11 . 
The  results  are 


[/ (£  -  *=l  -  /*?.)  fcFn(fcp) 

where  In  =  /*n  (/i?xr2  -  1  )Jf(l*ur9).  Let 


4-  c.c.  = 


2v  m0c2  d 
c/n  e  dt 


(*c  ±  iyc),  (9) 


(*.  ±*y„)  =  f  dk±i±(k±y<**’-“‘)  +c.c.. 


(10) 


We  substitute  Eq.  (10)  into  the  wave  equation,  Eq.  (9),  to  obtain  a  pair  of  equations  that 
relate  b^n  to  Keeping  in  mind  that  the  variable  of  integration  is  a  dummy  variable  in 
each  case,  we  solve  for  b+n  and  6_n: 

2i/  m0c 7  u>  —  v0k 


feTn(A:)  =  — i- 


:&(*)• 


(11) 


c/n  e  (o;2/c2  -  k2  -  p2u) 

Next,  we  need  to  solve  for  the  dependence  of  particle  motion  on  the  radiation  field. 
We  assume  that  x2  +  y2  <  r2  such  that  the  centroid  motion  depends  only  on  the  field  on 
axis.  The  beam  centroid  equation  of  motion  can  be  written  as 

^^2  ~  (*e  ±  iye)  +  ftqkqv0f{z)(xc  T  iyc) 

=  7*  +  V°jt )  ±  lo»)lr=°e_1U,t  +  C-C-I  »  (12) 

where  k\  =  2i//(/?272rJ),  and  fi,  =  eBq/j0m0c. 

Returning  to  the  quadrupole  field,  Eq.  (1),  let 


/(*)  =  5^/nCos(nfc9z). 


(13) 
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For  simplicity,  we  keep  only  the  n  =  1  term  of  the  expansion.  Substituting  (7)  and  (10) 
into  (12),  we  group  terms  with  the  same  electromagnetic  wave  number  k  and  evaluate  the 
fields  on  axis.  We  solve  for  in  terms  of  b±n. 


[-(«  -  vok)2  -  k]v20)  i±{k)  +  n,fct(/,/2K  [m*  +  kq)  +  u{k  -  jb,)] 


7o  m.0c 

We  substitute  (11)  into  (14),  to  obtain 

ftqkqfl 


C  € 

=  -i - z{u>  -  v0k)bTU(k). 


(14) 


where 


S{u,k)(±(k)  —  *_!  |fT(fc  +  kg)  +  (T(k  -  fc,)j  , 

5(„t)-(,K-^(-1  +  _ -JI.. -,  ) 


(15) 


(16) 


and  k%  —  2u/‘y0Iu.  Eliminating  we  find  that  the  longitudinal  waveguide  modes 
coupled: 

Tm(u ?,fc)£+(fc)  =  —  Sm-l£+(k  +  2kq)  +  Sm+  l£+{k  —  2fcg)j  , 

where 

Sm{v,k)  =  S(v,k  +  mkg),  (18) 


are 


(17) 


2m(w,fc)  —  Sm_15m5m+1  —  - K\  [Sm_i  +  Sm+1] , 


(19) 


K\  —  Kqkqfi/y/2  and  Kq  =  Uq/v„. 

From  Eq.  (17),  the  dispersion  relation  for  the  TEn  mode  in  the  presence  of  the  beam 
is  a  function  of  the  determinant  of  an  infinite  tri-diagonal  matrix  of  the  form 


det 


/*.. 

• . 

0 

0 

\ 

T-2 

%S-  3 

0 

^5+1 

T0 

Ki  q 

•  .  . 

0 

%s+3 

T+2 

V... 

0 

0 

**./ 

=  0. 


(20) 


It  can  be  shown  from  Eq.  (20)  that  the  growth  rate  of  the  instability  is  periodic  in 
k.  For  a  given  unstable  frequency  u>0,  the  unstable  wave  numbers  are  at  all  k  =  jfc0  +  nkq, 


6 


where  n  =  0, ±1,....  is  an  integer,  and  Jfc0  is  the  unstable  wave  number  associated  with  a 
vacuum  waveguide  mode.  Coupling  to  an  infinite  number  of  modes  may  be  avoided  in  the 
approximation  that 

4w<1-  <21> 

It  will  be  shown  below  that  2^*  <  1  is  required  for  phase  stability.  To  zeroth  order,  the 

Kq 

approximate  dispersion  relation  is 

To  =  0.  (22) 


To  first  order,  we  find 


/  r_2  £s_3  o  \ 

det  I  To  ^5_!  =0.  (23) 

\0  ^S+,  T+2 
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III.  Analytical  Beam  Line  Decomposition 


Much  insight  can  be  obtained  from  analytical  decomposition  of  the  approximate  dis¬ 
persion  relation,  Eq.  (22).  This  can  give  us  analytical  expressions  for  the  growth  rates 
and  instability  boundaries  in  parameter  space. 

The  dispersion  relation  can  be  rewritten  with  the  current  coupling  terms  grouped 


together  in  a  term  cr, 

W-iWW+iIIo  =  *0, 

(24) 

where 

Wm(u>,fc)  =  (w/c)2  -  {k  +  mfc,)2  -  /tjj. 

(25) 

nm(u;,  k)  =  ~  ^  [“m+l  > 

(26) 

am(u,k)  =  u)/va  -  (k  +  rnkq), 

(27) 

and 

<rm  =  klalal^al.AWrn+iW^  +  WmWm.x  +  WmWm^ 

~kl(Wm^  +  Wm  +  Wm+1)  +  kt} 

-\k\K'dWm  .  (28) 

The  polynomial  n0  can  be  rewritten  as 

n0(u;,fc)  =  (a0  -  *,)2A0  -  2kqK$a0,  (29) 


where 


Ao(“>,*0 


A  good  approximate  beam  mode  decomposition  can  be  obtained  by  dropping  the  last  term 
of  (29),  since  K%/k*  is  typically  much  less  than  1.  There  are  six  waveguide  modes,  i.e., 
W-\WqW+\,  and  there  are  six  beam  modes.  The  four  beam  modes  in  (30)  have  the  same 
characteristics  as  those  in  the  helical  quadrupole  case.2 

The  beam  line  decomposition  helps  to  define  the  approximate  boundaries  separating 
the  orbit  instability  and  the  electromagnetic  instability. 


8 


Orbit  Instability 

As  in  the  helical  quadrupole  case,  the  electron  beam  develops  orbit  instability,  when 

*J/4  <K2d,  (31) 

with  or  without  electromagnetic  modes  in  a  waveguide.  This  is  confirmed  by  the  numerical 
solution  of  the  complete  dispersion  relation  (23)  and  by  particle  simulations. 

Electromagnetic  Instability 

Electromagnetic  instability  exists  at  the  intersection  of  the 


W-x(  u,,*)  =  0 


(32a) 


waveguide  mode  and  the 

k  +  (-*,/ 2  +  y]k\! 4  -  K\)  -  u>/v0  =  0  (326) 


beam  mode.  This  instability  is  repeated  for  all  wave  numbers,  k  =  k0  ±  nkq,  where 
n  =  0,  ±1,  ±2, ...  When  the  intersection  does  not  exist,  we  find  that  the  beam  is  still 
unstable  for 


P  = 


(33) 


but  the  growth  rate  is  much  smaller. 

In  the  limit  of  zero  beam  current,  a  stability  diagram  based  on  Eqs.  (31)  and  (33)  is 
shown  in  Fig.  1 
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IV.  Analytical  Expressions  for  Growth  Rates  and  Group  Velocity 


Following  the  same  procedure  as  outlined  in  Ref.  2,  analytical  expressions  for  the 
spatial  or  temporal  growth  rates  of  the  electromagnetic  instability  can  be  obtained.  The 
derivation  for  the  spatial  growth  rate  is  presented  first.  Electromagnetic  instability  occurs 
at  a  frequency  and  wavenumber,  (u>a,ka),  which  satisfies  (32a)  and  (32b),  i.e.,  frequency 
ua  satisfies 

v/K/c)2  +kq  =  (wa/vo)  +  kq/2  -  yj k\ / 4  -  K\,  (34a) 

and  the  corresponding  wave  number  is 

K  =  \/K/c)J  -/*!i  +  *V  (346) 


The  derivation  of  the  spatial  growth  rate  is  given  below.  The  dispersion  relation  (24) 
can  be  rewritten  as 

(*!  -  *>)«*  +  k,f  -  *>)((*  -  k,f  -  k\) 

{( k  -  k2)\k  -  ktf  ~  -  k,)2  -  At|]  +  2k,Kl(k  -  fc,)}  =  (35) 

where,  at  u  —  u>a, 


».(fc)  ~kl(k  -  k,)\k  -  k,)2(k  -  k {(It2  -  kl)((k  +  k,f  -  kl) 

-kl{(k' -  kl)  +  ((k  +  k,f  -  kl)\]  (36) 

-  ^Ir-ik2  -  k2)((k  +  k,)2-  k2)(k  -  k,)2, 


fc0  =  y/{u>/c )2  -  /xfj,  Aj  =  u?/v0,  k2  =  (w/uo)  ~  kg,  k3  =  (u;/v0)  +  kq,  fc4  =  (u>/r0)  +  kq/ 2, 
AAn  =  y^l/4  -  Aj,  and  AA2  =  yJk\/A  4-  ifj.  Defining  fc  =  ko +kq+6k,  the  imaginary  part 
of  6k  is  the  spatial  growth  rate.  Requiring  fc<)  4-  kq  =  fc4  —  Ak\  and  using  the  approximate 
beam  modes  (32b),  the  zeroth  order  dispersion  relation  (35)  can  be  approximated  by 


pa6k2  4  pb6k  4  <t  =  0, 


(37a) 


where 


pa  -  4Afci(Afc|  -  Akl)k0(k  -  k2)2{k2  -  kl)((k  4  kq)2  -  fco)U=*.  .«■»=«.»  (376) 
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Pb  =  4kqK$(k  -  h)k0(k  -  fc2)2(fc2  -  fc2)((fc  +  kq)2  -  fco2)U=fc.,u,=u,.,  (37c) 


and 


&  —  ^0  |fc=fc.  ,d»=a>.  • 


(37d) 


The  derivation  of  (37a)  assumed  6k  <C  2fco  and  6k  -C  2Afci-  The  analytical  growth  rate 
expression  evaluated  at  (wa, ka)  is 


Im(6k)  =  y/(Pb/pa)2/4  -  <r/pa. 


(38) 


We  can  also  give  the  analytical  expression  for  the  temporal  growth  rates.  The  disper¬ 
sion  relation  (24)  can  be  rewritten  as 

^(a;2  -  u^jXw2  -  u>l)(u;2  -  a?2) 

j^(w  ~  ^4)2[(w  -  w,)2  -  Au;?][(u>  -u>qf  -  A w|]  -  2kqK$^(u>  -u>3)j  ~  /?®do,  (39) 
where,  at  fc  =  fc0, 

^o(w)  ^^|(o;  -  w2)2( uj  -  u>3)2(u>  -  w4)2  {  -^(a>2  -  w\){u2  -  u%) 
v0 

-  H[(*2_w?)  +  („2  _„*)]}  (40) 

-  ^  “  "’M"*  ~  w«)(w  ”  W2)2’ 


w-i  =  y/{k  -  fc,)2  +  pf1c,  u)0  =  y/k2  4-  p\tc,  u>i  =  y/(k  +  fc,)2  +  pfac,  u>2  =  {k  -  kq)v0, 
u>3  =  fcv0,  o>4  =  (fc  +  fc,)v0,  w,  =  (fc  -  kq/2)v0 ,  Awi  =  yJk2/4-  K\v o,  and  Aw2  = 

^/fc2/ 4  +  K%v0.  Defining  u;  =  u>_j  -f  6u>,  the  imaginary  part  of  6u?  is  the  temporal  growth 
rate.  Requiring  w_i  =  u>,  +  Au/j,  the  analytical  temporal  growth  rate  expression  can  be 
simplified  to 

Im  (^*9  =  V9^4  ~  9c 

where 

.  -  .  (ui  —  Ult  \ 

,  (416) 


W_i=W«,fc=fc, 


(41a) 


_  _  -405I,  if4  _ (w  ^3) _ 

Qb  P0«,Ad  AiVi(Au}2  _  A(jJ2^w  _  W4)2 


u>=u;_  1 
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and 


9c  = 


C10/3oO-  0 


4Au>j(Au;f  —  Au>i)o;_i(u>  —  o>4)2(u>2  —  a>o)(u>2  -  u)f ) 


(41c) 


u>=u;_  i 


From  our  numerical  solutions  of  the  higher  order  dispersion  relation  (23)  we  conjecture 
that  when  a  higher  order  dispersion  relation,  such  as  (23),  is  considered,  spatial  and 
temporal  growth  rates  at  other  intersections  in  (u>,  k)  space  of  waveguide  modes  and  other 
beam  modes  will  collapse  to  expressions  similar  to  Eqs.  (35)  and  (41).  (These  intersection 
occur  at  u>  =  u>0  and  k  =  ka  ±  ntcq). 

Because  the  beams  have  a  finite  length  and  because  the  group  velocity  of  the  instability 
is  slower  than  the  beam  velocity,  the  instability  can  propagate  out  of  the  tail  of  the  beam. 
From  (38)  and  (41a),  dropping  the  Pb  and  qb  terms,  the  group  velocity  of  the  instability  is 


vg  =  6u/6k 


“( 


“>«/c 


i/2y 


(42) 


The  number  of  e-folds,  N ,  that  can  occur  within  a  beam  of  length  tb  is 


(43) 


where  Le  =  1/Sk  is  the  e-folding  distance,  L  is  the  total  distance  traveled  by  the  beam 
and  Vb  is  the  beam  velocity. 
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V.  Numerical  Examples 


Here,  the  dispersion  relations  (22)  and  (23)  are  solved  numerically.  We  verified  that 

i)  Eq.  (22)  is  a  fair  approximation  to  (23), 

ii)  the  approximate  beam  mode  decomposition  Eq.  (29)  of  (22)  is  good, 

iii)  the  boundary  of  instability  is  as  predicted  by  Eqs.  (31)  and  (33)  and 

iv)  analytical  growth  rate  expressions  are  in  good  agreement  with  the  results  from 
(23). 

Parameters  for  numerical  examples  are  typical  of  a  high-current,  induction-accelerator 
beam,  It  =  1  kA  and  70  =  5.  We  chose  a  quadrupole  gradient  of  Bqkq  =  221  G/cm, 
/i  =  A/it,  quadrupole  wave  number  of  kq  =  0.5  cm-1  and  the  waveguide  radius  of  rg  = 
3  cm.  These  parameters  fall  in  the  electromagnetic  unstable  regime  in  Fig.  1,  where 
\/Vii/(7o  —  1)  =  0.125  cm-1  and  P  =  0.448  cm-1  (see.  Eq.  (33)).  Figure  2  plots  the 
(u>,  k )  diagram  corresponding  to  the  zeroth  order  approximation  to  the  dispersion  relation, 
Eq.  (22).  The  solid  circle  in  Fig.  2  marks  the  instability  denoted  by  Eqs.  (32a, b).  The 
instability  encircled  by  the  dashed  curve  is  an  erroneous  one.  It  is  modified  when  higher 
order  expressions  for  the  dispersion  relation  are  solved.  Figure  3  plots  the  first  order 
approximation  to  the  dispersion  relation  Eq.  (23). 

Figure  4  plots  of  the  spatial  growth  rate  as  a  function  of  wave  number,  k,  for  u  >  0 
from  (23).  The  instability  is  periodic  in  k.  The  growth  rates  of  all  the  instabilities  are 
identical,  except  one.  The  instabilities  at  (fc  ~  fc0  —  kq  and  k  ~  k0  —  2fcg)  are  from  coupling 
to  r+2,  where  k0  ~  0.22  cm  1 .  The  instabilities  from  coupling  to  T_ 2  (not  shown)  are 
at  (k  ~  ko  4-  2 kq  and  k  ~  ko  4-  3 kq).  The  instability  at  (k  ~  ko  —  2kq),  similar  to  the 
instability  marked  by  the  dashed  circle  in  Fig,  2,  is  erroneous,  and  is  modified  to  give 
the  same  growth  rate  as  that  at  ko  -f  kq  as  more  terms  of  the  full  determinant  (20)  are 
included. 

Figure  5  plots  the  spatial  growth  rates  as  a  function  of  frequency  u>/c.  for  u;  >  0 

for  i)  zeroth  order  approximation  from  Eq.  (22),  dashed  ( - )  curve,  ii)  first  order 

approximation  from  Eq.  (23),  solid  ( - )  curve,  iii)  exact  growth  rates  from  Floquet 

formulation,5  (+  symbols)  and  iv)  analytic  growth  rate  expression  from  Eq.  (38),  (• 
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symbols).  The  growth  rate  of  the  zeroth  order  approximation,  the  first  order  approximation 
and  the  Floquet  formulation  are  in  excellent  agreement. 

The  unstable  beam  mode  given  by  Eq.  (32b)  is  in  excellent  agreement  with  the  line 
in  (u>,  k)  space  in  the  dispersion  diagram  obtained  by  numerically  solving  the  dispersion 
relation.  Two  of  the  six  beam  modes  given  by  Eq.  (29)  are  not  as  accurate.  One  of 
these  two  beam  modes  appears  on  the  left-hand-side  and  one  of  these  two  beam  modes 
appears  on  the  right-hand-side  of  Fig.  2  Since  the  incorrect  beam  modes  are  not  the 
unstable  beam  mode  of  interest  and  are  sufficiently  far  from  the  unstable  beam  mode,  the 
analytical  growth  rate  expressions  are  in  good  agreement  with  the  numerically  obtained 
growth  rates  shown  in  Fig.  5. 

Figures  6,  7,  and  8  vary  the  quadrupole  gradients  while  keeping  all  other  parameters 

the  same  as  that  used  for  Fig.  2.  Figure  6  compares  the  numerically  ( - )  and  analytically 

(•)  obtained  spatial  growth  rates  for  three  different  quadrupole  gradients  (a)  Bqkqf\  =  200 
G/cm,  (b)  Bqkqfi  =  400  G/cm  and  (c)  Bqkqfi  =  600  G/cm.  Figure  7  compares  the 

numerically  ( - )  and  analytically  (•)  obtained  temporal  growth  rates.  The  comparison  of 

numerically  obtained  (hashed  area)  and  analytically  obtained  (•)  group  velocities  is  shown 
in  Fig.  8. 

The  boundaries  of  the  instability  regime  predicted  by  Eqs.  (31)  and  (33),  in  the  limit 
of  zero  beam  current,  are  in  good  agreement  with  numerically  obtained  boundaries. 

Numerical  verification  of  the  weak  unstable  regime  (V)  was  carried  out  by  varying  the 
energy  of  the  electron  beam  while  holding  the  rest  of  the  parameters  constant:  quadrupole 
gradient  of  Bqkqfi  =  200  G/cm,  quadrupole  wave  number  of  kq  =  0.5  cm-1 ,  the  waveguide 
radius  of  rg  =  1.5  cm  and  current  7e  =  1  kA.  For  7 „  >  2.85,  the  beam  is  in  the  unstable 
regime  (I),  while  for  y0  <  2.85,  the  beam  is  in  the  unstable  regime  (V).  The  boundary 
for  the  instability  is  given  by  Eq.  (33).  Small  growth  rates  are  observed  numerically  for 
7 0  <  2.85  from  the  first  order  dispersion  relation  (23),  (see  Fig.  9). 

For  larger  current  and  smaller  quadrupole  gradient,  the  difference  of  growth  rates 
between  regime  (I)  and  (V)  is  much  larger.  Figure  10  is  a  plot  of  the  temporal  growth 
rates  for  different  values  of  j0  for  the  parameters:  quadrupole  gradient  of  Bqkqf\  =100 
G/cm,  quadrupole  wave  number  of  kq  =  0.5  cm-1,  the  waveguide  radius  of  rg  =  1.5  cm 
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and  current  Je  =  10  kA. 

In  fact,  instabilities  at  the  intersection  of  unstable  beam  mode  (32b)  and  all  other 
waveguide  modes  also  exist  in  the  electromagnetic  unstable  region  (I),  except  that  their 
growth  rate  is  much  smaller,  (see  the  small  blip  at  a >/c  ~  0.672  of  Fig.  5). 
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VI.  Summary  and  Comments 


Electromagnetic  instability  on  an  intense  beam  in  a  FODO  lattice  was  studied  by 
deriving  and  analyzing  a  dispersion  relation,  which  is  a  matrix  equation  of  infinite  dimen¬ 
sions.  We  showed  that  numerical  solutions  of  an  appropriately  truncated  matrix  equation 
give  growth  rates  in  excellent  agreement  with  the  exact  growth  rates  as  determined  from 
Floquet  Theory.5.  Furthermore,  analytical  growth  rates  were  obtained  and  are  in  good 
agreement  with  numerical  solutions  of  the  dispersion  relation.  Stability  boundaries  in  pa¬ 
rameter  space  were  established  and  verified.  The  analysis  will  be  extended  to  consider 
cases  which  include  an  additional  nonzero  axial  (solenoidal)  field. 

The  nature  of  the  instability  is  such  that  growth  should  be  observable  for  high  current 
(>  100A),  moderate  energy  (7  <  100)  electron  beams.  Furthermore,  the  group  velocity  of 
the  instability  is  slower  than  the  beam  velocity  so  that,  in  general,  growth  will  not  be  ob¬ 
served  in  beam  pulses  of  less  than  several  nanoseconds  for  these  parameters.  For  instance, 
this  instability  is  not  knowhto  have  been  observed  in  rf  linear  accelerators,  presumably 
because  the  growth  rates  are  too  small  and/or  the  electron  beam  pulse  length  is  short. 

Ion  beams  are  not  expected  to  produce  electromagnetic  instabilities  because  ion  beams 
are  typically  very  slow,  i.e.,  (3  <  0.2.  The  parameters  of  typical  ion  beam  falls  in  the  weak 
unstable  regime  of  the  instability  diagram.  Damping  of  the  TEu  mode  could  easily  be 
accomplished  with  finite  Q  values  in  the  waveguide  or  the  accelerator  cavity. 

Any  high  current  induction  driven  electron  accelerator  utilizing  FODO  lattices  would 
be  susceptible  to  this  instability  under  the  idealized  conditions  stated  here.  Realistic 
effects,  such  as  energy  spread,  nonlinear  mode  mixing,  wave  saturation  may  significantly 
reduce  the  effect  of  the  instability. 
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Appendix  A:  Proton  or  Ion  Beams 


For  proton  or  ion  beams,  the  coefficients  of  the  dispersion  relation  are  modified.  We 
will  consider  ion  beams  with  velocity  Uj,  mass  m*,  charge  per  ion  g  and  ion  current  h  = 
geviTii(irrl),  where  n*  is  the  ion  density  particle  density  and  is  the  ion  beam  radius. 
The  complete  dispersion  relation  is  identical  to  Eq.  (20),  with  the  definitions  of  several 
constants  replaced  by  those  given  here: 


n,  =  —  b„ 

7771,-  C 


{A.  1) 


Kq  —  SI  q/vi) 

9 2e2  ,  geli 

v  =  - ~ni(irrb)  =  - 

miC 1  rriiViC * 


(A.2) 

(A.3) 


L  2  _  ^  l' 

Kfc  —  7 - • 

hi  7 

Dominant  instability  exists  if  the  beam  mode 


(AA) 


k  +  (-kq/ 2  4-  Jfcf/4  -KD-  u/vt  =  0 


(A.  5) 


intersects  the  W_i  =  0  waveguide  mode.  Weak  instabilities  exist  at  intersections  of  (A. 5) 
with  other  electromagnetic  waveguide  modes. 
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Quadrupole  Wave  Number  k  (cm 


ts|m121a7'2-1>  (cm_1 ) 

Fig.  1.  Instability  diagram  in  the  limit  of  zero  beam  current. 
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oj/c 


Wave  Number  k(cm-1) 

Fig-  2.  Dispersion  diagram  of  the  zeroth  order  approximation  of  the  dispersion  relation  Eq. 

(22). 
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(x)  /c 


-1.0  -0.5  0  0.5  l.o 

Wave  Number  kCcnrT1) 

Fig.  3.  Dispersion  diagram  of  the  first  order  approximation  of  the  dispersion  relation  Eq. 
(23). 
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patial  Growth  Rate 


0.64  0.66  0.68 

co/c  (cm-1) 


Fig.  5.  Spatial  growth  rate  as  a  function  of  frequency  w/c  for  i)  zeroth  order  approximation, 

dashed  ( - )  curve,  ii)  first  order  approximation,  solid  ( - )  curve,  iii)  exact  growth 

rates  from  Floquet  formulation,  (+  symbols)  and  iv)  analytic  growth  rate  expression, 
(•  symbols). 
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Spatial  Growth  Rate 
T  ( 1 0-2/ cm ) 


0)/c  (cm  ') 

Fig.  6.  Spatial  growth  rates  for  w  >  0  and  (a)  Bqkqfi  =  200  G/cm,  (b)  Bqkqfi  =  400  G/cm 
and  (c)  Bqkqf\  =  600  G/cm,  while  keeping  all  the  parameters  the  same  as  Fig.  2. 
The  analytically  calculated  growth  rates  are  indicated  by  (•). 
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Quadrupole  Gradient  Bq  kqf,  ( 100  G/cm) 

Fig.  7.  Comparison  of  numerically  ( - )  and  analytically  (•)  obtained  temporal  growth  rates 

for  the  same  parameters  as  Fig.  6. 


Group  Velocity  /c 


Quodrupole  Gradient  Bq  kqf,  ( 100  G/cm) 

Fig.  8.  The  comparison  of  numerically  obtained  (hashed  area)  and  analytically  obtained  (•) 
group  velocities  for  the  same  parameters  as  Fig.  6. 
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Temporal  Growth  Rate,  f/c  (1 


Fig.  10.  Temporal  growth  rate  versus  7,,  for  Bqkqf\  =100  G/cm,  kq  —  0.5  cm  1  r9  =  1.5  cm 
and  Je  =  10  kA.  The  instability  is  in  regime  (I)  for  7„  >  2.9  and  in  regime  (V)  for 
7o  <  2.9. 
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